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We analytically find an approximate Bloch-Messiah reduction of a noncollinear parametric am- 
plifier pumped with a focused monochromatic beam. We consider type I phase matching. The 
results are obtained using a perturbative expansion and scaled to a high gain regime. They allow 
a straightforward maximization of the signal gain and minimization of the parametric fluorescence 
noise. We find the fundamental mode of the amplifier, which is an elliptic Gaussian defining the 
optimal seed beam shape. We conclude that the output of the amplifier should be stripped of higher 
order modes, which are approximately Hermite-Gaussian beams. Alternatively, the pump waist can 
be adjusted such that the amount of noise produced in the higher order modes is minimized. 

PACS numbers: 42.65.Yj Optical parametric oscillators and amplifiers, 

42. 65. Re Ultrafast processes; optical pulse generation and pulse compression, 

42.65.Lm Parametric down conversion and production of entangled photons. 



1. INTRODUCTION 

Parametric down-conversion is an important and versatile source of light. Its applications span from am- 
plification of laser pulses to creation of strictly quantum states in the form of entangled photon pairs and 
squeezed vacuum. As optical parametric amplifiers rely on nonresonant nonlinear interaction, they can 
be utilized to amplify extremely wide bandwidth pulses. One of the most appealing applications of this 
phenomenon is the Optical Parametric Chirped Pulse Amplifier (OPCPA) |TJ |3] in which ultrashort pulses 
are stretched, amplified and then compressed again. The gain medium is typically pumped by a nanosec- 
ond second harmonic pulse from a Q-switched laser. Such amplifiers can provide a few femtosecond long 
multi-terawatt pulses [21 0]. 

One of the most important and fundamental problems with this type of amplifiers is the presence of spon- 
taneous parametric fluorescence, which unavoidably accompanies useful gain and compromises the pulse 
contrast ratio. This issue is central in many applications and has been extensively studied [5HZ]. Although 
the numerical results provided in [7] are very accurate, they are based on Langevin noise equations [H] de- 
veloped for a cavity parametric amplifier. Consequently, the noise parameter has to rely on an experimental 
calibration. Furthermore, their numerical character does not give direct guidelines for noise reduction. 

In this paper, we develop a simple analytical model of a 3D Noncollinear Optical Parametric Amplifier 
(NOPA). We assume type-I noncollinear phase matching in a second-order nonlinear crystal pumped with 
a monochromatic Gaussian beam. As a consequence, we are able to calculate the intensity of parametric 
fluorescence for any given signal gain, as well as its spatial and spectral distributions. Moreover, we find the 
optimal seed mode and pump waist for which spontaneous fluorescence is minimal. 

Our model is built around the Bloch-Messiah reduction [5] of an optical parametric amplifier operated 
without pump depletion. It is a theorem which allows us to find a set of orthogonal, characteristic modes (or 
eigenmodes) of the amplifier and their respective gains. In the case of monochromatic pumping, the modes 
are the input and output beam shapes of the signal and idler fields at any pair of frequencies matching up 
to the pump. The main feature of the reduction is that if the amplifier is seeded with a beam matching its 
characteristic input mode, then this beam is amplified by the associated gain factor and assumes the output 
characteristic shape on the other end of the amplifier. In addition, spontaneous parametric fluorescence is 
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produced in all of the characteristic output modes and the number of photons scattered into each of those 
modes is equal to the gain of that mode. In particular, the immediate conclusion from the reduction is that 
the optimal signal to parametric noise ratio is obtained when the seed beam has the shape of the fundamental 
input mode. The output of the amplifier may be spatially filtered to reject as many of the higher other modes 
as possible. 

The paper is organized as follows. Section 2 presents general theory of parametric down-conversion, with 
emphasis on the 3D crystal pumped by a monochromatic wave. We introduce the Bloch-Messiah reduction 
theorem and discuss its general consequences. In Section 3 we derivate the characteristic modes and their 
gains for two particular settings - noncollinear with no walk-off and noncollinear with a significant walk-off. 
Section 4 summarizes the results of the analytical reduction. In Section 5 we compare analytical results with 
precise numerical calculations and discuss the accuracy of the approximations assumed. Section 6 describes 
behavior of the parametric gain in the intense pumping regime, provides estimation of the fluorescence noise 
in practical situations and gives solutions for the noise reduction. Finally, Section 7 concludes the paper and 
gives insight into possible applications and extensions. 



2. BLOCH-MESSIAH REDUCTION 



We aim at investigating the properties of a 3D parametric amplifier pumped by a strong, monochromatic 
laser beam of frequency uj p and diameter 2w p . We assume that the signal and the spontaneous fluorescence 
fields are always so weak that they cannot influence the evolution of the pump. Within this approximation, 
the nonlinear interaction with the pump field couples pairs of conjugate frequencies, commonly called signal 
uj s and idler = u> p — u> s . Let us note that any two adjacent frequencies w s ^ lo' s evolve completely 
independently thanks to narrowband pumping. Therefore, amplification of a broadband seed pulse can be 
calculated separately for each of its monochromatic components, which may enter the crystal at different 
angles [lOj . Throughout the paper, we will analyze the interaction of a particular pair of frequencies w s and 
ixi i with the pump. 

We choose a coordinate system where the z axis lies along the central fc-vector of the pump k p and the crystal 
faces are perpendicular to this axis. Then, it is most convenient to follow the convention of describing the 
evolution of the system as a function of z. We expand the electric field of each of the interacting waves in 
the basis of monochromatic plane waves parameterized by frequency u> and perpendicular components of the 
wave vector k x and k v 



E(t,x,y,z)=i du;dk x dk y . ^^ e -*^+***-+*if*i>A(a; J k x , k y , z) + ex., (1) 

J y 47T£ cv n ( a; ) 

where A(lo, k x , k y , z) is the z-dependent frequency-space amplitude of the wave. With the above normaliza- 
tion k x , k y , z)\ 2 is the photon density. 



2.1. 1-D amplifier 



Before analyzing a bulk amplifier let us briefly recall the basic concepts for a waveguide amplifier which 
can be reduced into a 2-mode amplifier for a defined signal frequency. This simple system will help in the 
comprehension of a full-fledged 3-D amplifier. Later on, the Bloch-Messiah decomposition will allow us to 
simplify the 3-D amplifier to a set of independent 2-modc amplifiers. 

The pump field undergoes only dispersive evolution and its amplitude along the z axis is A p (z) — Pe lfcp zZ , 
where k pz is the pump wavenumber and P is the initial amplitude. The evolution of the signal amplitude 
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A s (z) and the idler amplitude Ai(z) is described by the equations of motion 

~A s (z) = ik s , z A s (z) + X Pe ik ™ z A*(z) 

~Mz) = ihzAi(z) + x Pe lk ^ z A* s (z). (2) 

The wave vectors for the signal and the idler are k SiZ and ki tZ , respectively. The first term on the right 
hand side is responsible for linear propagation. The second one describes the nonlinear interaction with the 
coupling constant x- 

For a waveguide which starts at —L/2 and ends at L/2 one can obtain the following input-output relations: 

(e l53 A° ut ) = cosh(()(e* 5l Al n ) + sinh(£)(e ll5 Mf)* 

(e l5 *A° ut ) = cosh(0(e J,52 An + sinh(£)(e l51 A\ n )* . (3) 

The upper indices in and out denote positions z = —L/2 and z = L/2, respectively. Phase terms with Si, 
82, S3 and ^4 depend on wave vector matching and we skip their explicit form for brevity. The intensity gain 
of the amplifier equals cosh 2 £, where the gain parameter £ is proportional to the pump amplitude P. In 
addition to the amplified signal, the interaction produces spontaneous fluorescence. Quantum-mechanical 
considerations [H] show that the average number of photons generated this way equals sinh 2 £. This system 
is called a 2-modc or nondegenerate parametric amplifier [13 . If one needs to calculate the amplification 
of a realistic laser pulse, the pulse has to be decomposed into monochromatic waves via Fourier transform, 
amplified as in ^ and composed back. 



2.2. Bulk amplifier 



Let us now switch to a bulk crystal with the length L in the z direction and infinite in the x and y directions. 
The pump A p (k p ^±_), the signal ^4 s (fc s and the idler A^fc^jJ amplitudes have spatial freedom which is 
encoded in their dependence on the spatial wave vector k± = (k x , k y ). The infinite transversal size of the 
crystal allows interaction only of those signal and idler components for which transverse wave vectors sum 
up to the pump's. Therefore in each slice of the crystal with given z only those interactions can take place, 
which preserve both the energy lj p — uj s + Ui and the perpendicular components of the momentum 

k p ,± = fc s ,_L + fcj.x- (4) 
Using these principles, an equation of motion for the waveguide amplifier |2| can be extended as follows 

g z A s (k St± ,z) = ik SiZ A s [k s , x ,z)+X / dk h± A p (k s ^ + k itl _)e^-' z A*{k h± ,z) 

^A i (hA-,z) = ik i , z A i (k i>A _,z)+x J dk s ^A p {k s . x + ki^)e lk ^ z A* s {k s ^,z). (5) 

As in the waveguide |2]), the first term on the right hand side represent the dispersive propagation. The 
second term is responsible for nonlinear interactions involving all the possible planewave components fulfilling 
Q. The above equations are linear in the signal and idler amplitudes, and so are the input-output relations 
for the amplifier. Formally, they can be written down using input-output relations with integral kernels (or 

Green functions) C ss (fc s ,j_, fes.x), S a i(k Sy ±, Cu(k it j_, ki,x), S is (k it x,ks,x) 



A° s ut (k s>1 _) = J dk l s ^C ss {k s ,x,K^)AT{K.x) + J dk i>± S si (k St± ,k i>± )A} n *(k i>± ) 

A° ut (h±) = I dk^Cuik^k'^AY 1 ^) + f dk 8 , ± S is (ki, x ,k 8>x )Af*(k i>± ). (6) 
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Since the integral kernels describe a reversible process, they turn out to be interdependent. The Bloch- 
Messiah theorem [5] relates their singular value decompositions (SVD): they have common modes and gain 
parameters, that is 

oo oo 
C S8 (k s , ± , k' s ,x) = E C*&,X) COBh(CnX n *(^,x). C tl (h±, ki, X ) = £ Cfox) COsh(£ n )C*(^ x ), 

oo oo 

- C"*(^,x) sinli(en)C(^,x), 5 is (fc i>± , fc s ,x) = E sinh(e„X"(fc s ,x). 

n=0 n=0 

(7) 

Above, ip™(k Sj ±), VC'(^s,_l)j and (**,!.) represent the input signal modes, the output signal 

modes, the input idler modes and the output idler modes, respectively, while £ ra are their gain parameters. 
Mathematically, ^™(fc S) j_), tf'^Qis ±), </>™(fci,_i_) and ^"'(fcj.j.) are four different sets of orthonormal func- 
tions. Physically, they represent characteristic beam shapes in the far field. Modes with different n indices 
do not couple. More precisely, the input-output relations given by ^ can be brought to the canonical 
Bloch- Messiah form if we decompose the signal and the idler amplitudes in the bases of their characteristic 
modes 

oo oo 
oo oo 

Ar(h,x) = T / ^n(kx), Ar t (k i ,±) = J2 A ^>n t ikx)- (8) 



n=0 n=0 



Substituting the above into (|6j) and using Q we find that the whole OPA can be decomposed into a set of 
independent 2-mode amplifiers 



= cosh(£ n )A™ + sinh(C„)^ 
Aft = cosh(£„) A% + sinh(£„)A™* , (9) 

where n indexes subsequent independent modes, conventionally sorted by their gains. Note that the above 
equations have the same form as the input-output relation for the waveguide 2-mode amplifier given in ([3|. 
Similarly, the intensity gain of the n-th mode equals cosh 2 £„ and the average number of the photons spon- 
taneously scattered into that mode equals sinh 2 £„. Therefore, we are interested in finding the characteristic 
beam shapes ip™(k s< ±), ip^(k a ±) and gain parameters £„, as this will provide the complete information. 
In particular, seeding the amplifier with its fundamental mode V>q"(A; S) _l) would give us an advantage of 
obtaining maximal possible gain cosh 2 £o- The idea of the decomposition is illustrated in Fig. [l] 

Let us make a nontrivial observation. For a special case of the crystal surfaces equally separated from the 
z = plane and with the real pump amplitude in z = 0, we get a symmetry S S i(k St ±, ki t ±_) = S* s (ki t _L, k St ±) 
[13] • When every gain parameter is unique ^ £ m , the input and output signal modes are mirror reflections 
of each other, as are the input and output idler modes 

This observation significantly simplifies the Bloch-Messiah reduction given in Q. In particular it allows us 
to deduce all the characteristic mode functions from the SVD of Si S (ki y ±, k s ,±)- 



Even though in and out modes are similar (10 1, the behavior inside the crystal may be complex. Our 
description says nothing about what happens between the surfaces - we have only the input-output relations. 
The shape of the modes and their amplification vary in a number of parameters, including the crystal length 
L. 
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FIG. 1: Schematic of a NOPA pumped with a focused monochromatic beam. When an incident signal beam is in the 
i/)™(fc 3 ,x) mode, it produces the signal (fc s ,x) at the output with the amplitude amplification of cosh(£ n ) and 
the idler (j)™* (ki t ±) with the amplitude sinh(£ n ). 



2.3. Perturbative expansion 



In general the integral kernels defined in (|6| cannot be found analytically. However, one can obtain closed 
form expressions in a low gain regime x — > 0. This is accomplished in two steps. First we assume zero 
nonlinearity \ = 0. The ([5| can be immediately integrated as follows 

A s (k St x,z) = A s (k S:± ,z = 0)exp[ik s ,z(k B ,±)z], Ai(k it x,z) = A t (k ii± ,z = 0) exp[ik itZ (k it x)z], (11) 

which describes the linear propagation. To obtain the first-order approximation with respect to x, we need 
to integrate ^ over —L/2 < z < L/2 with the zero-order approximation (111 substituted on the right hand 
side. Comparing the result to the definition of the integral kernels (JsJ) , we obtain 



C ss (k s ,x, k s ,x) = ex.p(ik s>z z)S(k' s z - k s<z ), 



Cu(ki,±,ki t ±_) = exp(iki >s z)S(k' i , - h. z ). (12) 



In the first-order approximation the above functions describe only linear propagation. However, the remain- 
ing two kernels 



S si (k s ,±,ki,±) = X Le~^~ k °^ L/2 A p (k hl _ + fe s>± )sinc (AfeL/2) 
S u (ki,±, fc.,_i_) = X Le t{k ^- k ^ )L/2 A p (k Zi± + fc s , ± )sinc (AfcL/2) , 



(13) 



describe down- conversion and give us insight into the characteristic modes of the Bloch-Messiah reduction 
([7]). We call Si s the scattering kernel, as it describes amplitudes with which the incident signal creates the 
'scattered' idler beam. The Afc is the wave vector mismatch along the z axis in this process 

Afc = k s , z + k hZ - k p , z . (14) 

The z components of the wave vectors are functions of their transverse wave vector k SyZ = fc s ,z(fc s ,_L), 
ki t z — ki^ z (ki x) and k pz = kp Z (kix + k s j_). 



Once we perform SVD of the scattering kernel Si S (131 we find all the characteristic modes and gain pa- 
rameters of the Bloch-Messiah reduction ^ with the help of the identity (10). Obtaining an approximate 
analytical SVD of the scattering kernel SZ is described in the next section. It requires several nontrivial 
steps. However, we will always separate the variables in cylindrical coordinates as follows 



Sis{h,_L>ipi\ fc s ,j_, tp s ) = K(k it x,k s ,x)®(<Pi, <Ps) 



s,_U 



; I ; -in 



n m " »- *N-v~*s . ' 



(15) 



G 



where functions and singular values with upper indices K and <£• are singular value decompositions of 
K(h,±,ks,±) and $>(<pi,(p s ), respectively. Above we express the SVD of Sj s as a product of SVDs of its 
radial and angular parts, K(ki±, k Sj ±) and <&(cpi,tp s ). Note that this is also an approximation which relies 
on the fact that the radial mode functions a^ OMt (fci.j_) and ipn' m (ks,±) are centered around certain k± ~ k 
and do not extend towards the origin kj_ ~ 0. 



3. ANALYTICAL APPROXIMATION 



Our task is to find the SVD of the scattering kernel Si S (ki y ±, h s ,±) given in (131 for a Gaussian pump beam. 
We will accurately find a mode with the highest amplification. To obtain the result we go through several 
approximations as illustrated in Fig. [2j They include replacing sine with a Gaussian, Taylor expansion of 
the phase mismatch AkL/2 and the separation of variables. We will consider two different physical settings 
- noncollinear without pump walk-off and noncollincar with a significant pump walk-off. For each of them, 
we obtain closed form Bloch-Messiah reduction, as summarized in Section 4. We restrict ourselves to type I 
down-conversion, although we conjecture that similar approach may work also for type II process. 
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General theory 
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FIG. 2: Schematic of the procedures applied in Section 3 along with parts of Section 2 and Section 4. Solid lines 
stand for equalities, while dashed lines for approximations. Final results are summarized and discussed in Section 4. 

We choose Gaussian pump beam profile with waist w p and the wavefront parallel to the crystal surface. The 
pump amplitude in spatial frequency domain is 



Ap{kp^xi kp.y) — P i— exp 
where P is an amplitude factor. 



k 2 x 

p,x ' "p^yi 



(16) 



In the first step we approximate sinc(A/c/L) with a Gaussian. This erases the oscillating tails but enables 
further calculations ITS]. We write 



sine 







' (AkL) 2 ' 


(T) 


~ exp 


20 



(17) 
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bearing in mind we chose arbitrarily the width of Gaussian fit as illustrated in Fig. [3j Note that the length 
of the crystal L appears only in the argument of sine in ( 13 ). Consequently, the uncertainty in the numerical 



factor in the above approximation can be understood as a ~ 10% uncertainty of the crystal length L. 
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FIG. 3: Plot of sinc(a;) (solid line) and exp(— x 2 /5) (dashed line). In the approximation of sine with exp(— x 2 /k), 
reasonable values of the Gaussian width coefficient k spread from 4.16 to 6.28. One may use the closest Gaussian in 
the norm distance (k = 4.16), set the same FWHM (k = 5.18), compare Taylor series up to the second order (k = 6) 
or guarantee the same volume (k = 6.28). In various papers [14H18] different k values were used. We arbitrarily 
adopt ft = 5. 

Further steps involve Gaussian approximation of the scattering kernel. The task is simple but not straight- 
forward. We need to choose the suitable variables as well as the right Taylor series approximation. 
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FIG. 4: The crystal in a coordinate system. In both x and y direction the crystal is infinite, while in z its length is 
L. The optical axis lies in the xz plane, at the angle 9 to the z axis. 

Since the amplifier discussed has nearly cylindrical symmetry, it is convenient to introduce cylindrical coor- 
dinates in the following form 



[k a ± cos(ip s ),k s>± sia(<p a )] 



-ki t cos(ipi), 



sm(tpi)} . 



(18) 



Note that angles on the parametric down-conversion cone <p s and (fi are measured from 0° and 180°, re- 
spectively (Fig. E). Then the phase matching is achieved for tp a sa tpi. In cylindrical coordinates with 
approximation (17 1 the scattering kernel (131 reads 



S is (k itX , k s , x ,<Ps) = x PL^e^- k ^ L / 2 exp 

V 71 " 



fc? , - 2k s ,xki,± COs((f s - (fi) 1 A , 2 2 



2/^ 



20 



-Ak'L' 



(19) 



The scattering kernel achieves its maximal absolute value when k St ± + ki t ± — and there is no mismatch 
Afc = 0. We will linearize Ak around this point, which physically corresponds to the cone of ideal phase 
matching. 



To calculate the phase mismatch Ak given in ( 14 ) we need expressions for the z components of the wave 



vectors of the interacting waves. We consider type I down-conversion with the signal and the idler propagating 
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as ordinary waves and the pump as extraordinary. The optical axis lies in the xz plane, at the angle 9 to 
the z axis. We find that 



w 2 nl(uj s ) 



-k 2 . 

s,_L 



(20) 



(uj p - uj s ) 2 nl(uj p - uj s ) 



- k 2 



( k s,x "h k i,x) 



sin 



(29) [n*(w p )-nl(u p )] 



n (ui p )n e (oj p ) 



2 [nl(w p ) sin 2 ((9) + n 2 (uj p ) cos 2 {9)] n 2 (uj p ) sm 2 (9) + n 2 (uj p ) cos 2 (9) 



^ ( k s,y "H ki,yj 
72" 



[n 2 (u p ) sm 2 (9) + n 2 e (w p ) cos 2 ((9)] - (k s . x + fc^) 2 



where c is the speed of light, n (oj) and n e (uS) are ordinary and extraordinary refractive indices, respectively. 
The lengths of the wave vectors of the interacting waves are k s o, fc^o and k p and p is the tangent of the 
pump walk-off as illustrated in Fig. [5] 



k s fi = —n (uj s ) 



c 



«p,0 



n o (i0 p )n e (uj p ) 



n 2 {u p ) sin 2 (0) + n 2 (u p ) cos 2 (6») 



P 



- n (uj p - uj s ) 

sin(26>) [n 2 (w p ) - n 2 e (uj p )] 
2 [n 2 {uj p ) sin 2 [9) + n 2 e (uj p ) cos 2 (9) 



(21) 




FIG. 5: Schematic of ideal phase matching. The signal and the idler interact with the pump central component. 
Their respective wavenumbers are fc s ,o, fci,o and k Pt o. The perpendicular wave vector component of both the signal 
and the idler is fco. Inside the crystal, the signal travels at the angle arctan(a s ) w a s and the idler at arctan(ai) w a*. 

The phase mismatch vanishes Ak = when tp s — (pi and the perpendicular wave vectors components of the 
signal and the idler are equal k s ,j_ = fcj j_ = fco with 

^° = 2k o V^^°^° ^p-°^'° ^o^o) ~ ^io ~ — ^p,o- (22) 

The linearization of the Afc given in (14) around its zero with respect to fc Si j_, fci,j_ an d (</? s — </?,:) using (20) 
gives us 

Ak a [o, - pcos(^±^)] (fc fl , ± - fc ) + [a, + pcost^^)] (k %1 _ - k ) - k psm(^p±)(<p s - <p t ). (23) 



Both a s — fco/ k 2 — k 2 and a, = fc'o / \J k 2 — fc§ have an easy interpretation. Each of them is a tangent of 

the angle at which the signal or the idler travels inside the crystal in an ideal phase matching setting. Since 
p (21) is the tangent of the pump walk-off angle, the sum and difference of a s or at and p reflect how the 
or the idler diverge form the pump. 



signa 



The Gaussian pump (16) may be approximated by Gaussian function also in cylindrical coordinates (18). 
For the exponent we obtain 

k l± + k l± - 2k s ,_Lki,± cos(ip s - <pi) w k%(ip s - if,) 2 + (fc s ,_L - k it ±f, (24) 
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which is valid as long as l/w p <C ko, that is the down-conversion cone opening angle is much bigger than 
the pump beam divergence. 



Consequently, the scattering kernel Si s given in (19 1 substituted with (23) and (24 1 becomes a Gaussian 
function in k Sl ±, k^± and (tp s — ip^), that is 



P_ p i(k iiZ -k s ,z)L/2 



exp (—x Ax) , 



(25) 



where x = [(ip s — ifi), (k s ,i_ — ko), (fci,_i_ — fco)] is the deviation from ideal phase matching. The quadratic 
form coefficients are 



A = 



A A A - 

s*-ipip -f^-ips s^-tpi 

A A A ■ 

s^-sip -fi-ss ^-sl 

A- A A- 

^-lif ^HS 



^ + ^fc Vsin 2 (^) 








(26) 



the angular part 



T^ + f^K -pcos(^)] 2 
+ [a. - pcos{^p±)} [a, + /9C os(^4^)] 



+ Iff [a- ~ Pcos(^)} [a, + pcos(^)] 
^ + M^ + PCOs(^±^)] 2 



20 



fcopsin(^i) 





-[a s -pcos(^±^)] 
- [at +pcos(^4^-)] 



the radial part 

-[a s -pcos(^4^ 



- [a, +pcos(^4^)] 





the mixed part 



The matrix A is presented as a sum of different parts to facilitate the separation of variables. We call those 
components the angular part, the radial part and the mixed part. Once the mixed part is zero, the scattering 
kernel Si S may be written as a product of two kernels, as outlined in (15). We will develop two different 



approximations: one for no pump walk-off p — and the other one for a significant pump walk-off. The aim 
of further steps is to neglect the mixed terms and get rid of trigonometric functions. 



3.1. No walk-off 



In some special situations there may be no walk-off (p = 0). This is particularly easy to solve, as the quadratic 
form matrix ( 25 ) can be factorized into the angular and the radial part (as the mixed part vanishes) 



A = 



k O w p 

2 








2 ' 20 s 
w p L 2 
~2~ ' 20 a s a i 







20 
L 2 



a s ai 

„2 



20 "t 



(27) 



As the scattering kernel may be factorized Si S — K(ki t ±, k Si ±)$>(ipi, (p s ) as in ( 15 ), we need to find the SVDs 
of two separate parts. To deal with the angular part we write 



$(¥>»>¥'*) = CX P 



l (<Ps ~ ViY 



1 1 

/2~7r k w p 

1 1 

/2~7T k w p 



dm exp 



m 



2kW +im(iPs ~ iPi) 



exp(— imipi) exp 



m— — oo 



exp(-Mm<^ s ). (28) 
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We applied the Fourier representation of a Gaussian function and then we approximated the integration 
by summation. The second step required a moderately slowly changing function, which is guaranteed by 
l/w p -C ko, already required in the previous step (24 1. 

The remaining radial part of the kernel K(k{ t ±, k s> ±) is a quadratic form of (fc s ,j_ — fco) and (&i,_L — ko). It 
can be directly decomposed with the Mehler's Hermite polynomial formula [19 1 . which reads 



exp 



2(1 - M 2 ) 



(z 2 + y 2 ) + 



2fixy 



= y^V 1 - fi 2 fi n u n (x)u n (y), 



71=0 



u n {x) = (2"n!)-3ff n (x)exp(-x 2 /2), 



(29) 



where x, y and p are real numbers, while u n {x) are the Hermite-Gaussian modes. Comparing the coefficient 
of the radial part K(ki_j_, fc s ,j_) found in the lower-right part of (27) with the general Mehler's formula, we 
get 



K(k itX ,k s , x ) = xPL^e 



exp [-A ss (k s> ± - k ) 2 - 2A si (k St ± - k )(k it ± - k ) - A H {k itl _ - k ) 2 ] 



E 



x < e 



J c »fc 3 ,.L/2 u " [WiAki^J- ~ kp)} 
:k 3 ,zL/2 u n [w s ,r(k s ,A_ ~ kg)} \ 



2 n n 



x pl p yT^//-/< 



(30) 



where p, is a singular value scaling factor while w s and are width parameters for the signal and the idler, 
respectively 



\^~Ass~Aii / A ss Ai 



I A, 



ft', 



, P = y^2A7 s ( 1 



A 2 ^ 1/4 



1 



A A- 



,r — y/lAn ^1 



-4 2 



4 4 



1/4 



(31) 



The coefficients of the quadratic form matrix A are to be taken from (27). As outlined in (15), once we have 
the SVDs of both K(ki^±, k s j_) in (30) and 3?(<£j, f s ) in (28), we get the resulting Si S decomposition for the 
no walk-off p = case 



oo oo 



en, , \ ST ) iki ,L/2U n [wi, r (ki i± - k )} exp(-irmpi) 

V *-/Wi r 



n— m— — oo 



Dr 72^ v /l-/i 2 „ 
x xPL- /« exp 



rr; 



2fc>2 



-ik„ iZ L/ 



27T 



i/2 M„ [w s , r (fcs,_L - feo)] exp(OTl<£ s 



2tt 



(32) 



The resulting modes in a no walk-off setting are optical vortices. We discuss this result in Section 4. 



3.2. Significant walk-off 



The next case of our interest is the setting with a significant walk-off, i.e. when the total walk-off at the end 
of the crystal is of the order of the pump diameter Lp w w p . Again, we take the scattering kernel Si S given 
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in (25). The quadratic form A has nonzero mixed part, which prevents us from the separation of variables. 
In addition, it contains trigonometric functions of the sum of angles p s + pi. 

Let us first consider the physical picture of the significant walk-off setting. An efficient amplification of a 
beam may be obtained when it travels along the pump. It happens when tp a + (pi = or p s + (pi = 2n which 
corresponds to either signal or idler propagating along the pump. As the distinction between the signal and 
the idler is illusory, we may assume that the signal is propagating along the pump's walk-off, approximating 
all expressions containing (p s + (pi around 0. We will hold the quadratic term in A vtp and only the constant 
term everywhere else in A. Consequently, terms A vs — A atp and A V i = A^ disappear. Thus the quadratic 
form reads 



k o w p . Lll. 2 n 2 (* 
2 + 20 ^oP \~ 











20 y Us 



p? 



foK -p)( a i + p) 



o 



i( a » - p)( a * + p) 



20 \ a i 



Pf 



(33) 



Again, we have separated the variables Si s = K(ki.±,k s ^±)<fr(pi,p> s ) according to the scheme given in (15) 



The radial part K(ki.±, & s ,jJ decomposes with the Mehler's formula (29) as it did for the no-walk-off case 



All the coefficients are the same as those given in ([31]), with the quadratic form coefficient matrix A taken 
from a significant walk-off approximation matrix p3L 

The angular part <&{pi,p s ) contains not only (p s + (pi) 2 but also (ip s + (pi) 2 ((p s — ipi) 2 terms. They are 
difficult to handle, so we try the following approximation, with real variables p, q, x, and y 



exp 



p 2 (l + q 2 y 2 )x 2 ] w exp (-p 2 x 2 



p 2 q 2 (x 2 )y 2 ) = exp 



2 2 1 2 2 

-P x - -q y 



(34) 



That is, instead of the x 2 y 2 term, we took y 2 times mean x 2 , averaged over the exp f— p 2 x 2 ^j distribution. The 
rough approximation ( 34 ) is numerically checked to produce similar characteristic modes and gain parameters 
as long as q 2 <C 1. In our case q 2 = L 2 p 2 / (AOw 2 ) , so the total walk-off cannot be much larger than the pump 



width Lp <C V AOwp. With the approximation (34J) the angular part $(pi, p> s ) is a quadratic form of (p s and 
(pi. Utilizing once again Mehler's formula (29) we decompose $(pi,ip s ) in the basis of Hermite- Gaussian 
modes 



$((fi,ip s ) = exp 



h 2 ll\ 2 T 2 n 

K W 2 L p 
Tr^Ws-tPi) - 



Ps + Pi 



oo 

£ 

m=0 



u m {k w,p(pi 



k w„ 



(koWpips) 
y/l/(k w v ) 



(35) 



where p, v is the singular value scaling factor while fc u> v is the angular width parameter for both the signal 
and the idler 



(i lp =[l- 



1 



Lp 



2%/TO k w 2 



I 1 + 



1 



Lp 



2%/lfJ k wl 



(36) 



Hence, we obtain the SVD of the scattering kernel Si S in the significant walk-off approximation (33) with 
separated variables (fl5|), which reads 



oo oo 



Si s {k it ±,(pi;k s ,_L,(p s ) = 2 / 2 / ) c 



ikj.zL/2 U n [Wj,r(ki,± - hp)] U m (k W v (pi) 
y/l/w^r v /l/(fc ^ v ,) 



w. 



x X PL- 



' 2 'I-Pl 



koWpy/W^rWi 



-« e 



-ik.^L/2 U n [w s ,r{ k s.X ~ fc o)] u m(k W v p) s ) 
\/l/Ws,r \A/( fc C)"V) 



(37) 



Further discussion of the modes is given in Section 4. 
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4. CHARACTERISTIC MODES OF THE AMPLIFIER 



In Section 3, we applied a series ol approximations to obtain the singular value decomposition of the scatter- 
ing kernel Si s (13). The SVDs for the no walk-off ([32]) and significant walk-off (37) settings were compared 
with the general Bloch-Messiah reduction We identified the signal f/;° u 4(fc S; _L, 



the idler ^(fcj.x, ipt) = #£,(k 



><Ps) = tp l ^{k s ,±,(p s ) and 
_, <fi) modes - see |7| and Fig. [I] In particular, the highest achievable am- 



plitude amplification £ ,o is obtained with the seed tpoo{k a ,±,<Ps) on the input, which produces ^o^o (k s ,±, VsJ 
on the output. It needs to be remembered that the amplitudes in the spatial frequency domain are effectively 
the far-field image. 



4.1. No walk-off 



In the setting with no walk-off p = the modes resulting from the solution ( 32 ) are optical vortices 

zb out *(k , -a,)-**" (k , x ) - c-'*"^ 2 "" Klfe - fc ")] exp(im^ s ) 

i/l/W Sl r V2ir 

jjmt*(U „\ Ain (u ,„\ -ikj »L/2 M - K,r(fci,± ~ M] exp(-ZTO^) 

y/l/Wir V27T 



The minus sign at the </j s and (pi angles is a nontrivial consequence of twofold degeneracy of the singular 
values, which changes relations (10). The above modes are trigonometric functions on the down-conversion 
cone as plotted in Fig. [6] Their shape in the radial direction is that of a Hermite-Gaussian function with 
the width parameter 



W s ,r = y/l^V ^V w p + ^Q a s w i,r = y/ 1 - v 1 \j w p + ^Q a h ( 39 ) 

for the signal and the idler respectively, where v — {w^ + fo a s) ( w p + To a ?) / ( — w p + io a s ai ) ■ The gain 
parameters are 



|2 



j n ,n = XPL v x M p n exp ( - ) , = 07-7^7^1. (40) 

The characteristic modes are optical vortices 20J, thanks to the exp(zm</J s ) factor. Each of them carry orbital 
angular momentum, hm per photon. As the spontaneous down-conversion creates entangled photon pairs 
with the opposite angular momenta |21j , they may be used in the field of quantum cryptography |22j . The 
advantages of using entangled vortex states as media for quantum information include easy measurement 
[2"5] and the possibility to create Hilbert space of an arbitrary dimension [53] . 

To obtain the above decomposition we assumed that: l/w p <C fco> l/^s,r ^ &o and l/wi ir "C &o- This is 
approximately equivalent to the requirement that beam diffraction angles are much smaller than the angle 
between signal or idler and the pump. 
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FIG. 6: Three different modes for the no walk-off case: ipofi(ks,±,fs), ipo%(k s ,±, <fis) and ipl™4(k s ,±, <p s ). The phase 
e -»fcs,zi/2 ' g no ^. snown "jjjg figure represents the real part of the modes which, due to degeneracy, are also proper 
characteristic modes. 



4.2. Significant walk-off 



In the setting with a significant walk-off p > the modes resulting from (37) are elliptic Hermite-Gaussian 
beams 

lh out *(k , in\-lh in (k , CO ) - c ^k s ,zL/2 U n [w s , r {k s ,± - fcp)] M m (fc W y y s ) 
Vn,m l K s,_L, Ps) — V n .m\ K s.X^s) — C , , r 

VVV \/l/( fe o«v) 



±out*(u ,~\ J.in /;„ ,„ \ „-iki ,L/2 U n[ w i,r{k 



»,± - fco)] u m (k w v ipi) 



(41) 



as plotted in Fig. [7] That is, the characteristic modes are 2D Hermite-Gaussian functions, curved on the 
down-conversion cone. Their width parameters are 



= 0^=1^3 + ^(o a -p)a 
1 



£2 

i, r = \/l - V^XjUll + — (dj + p) 2 



ao 



V Lk oP 



[-wj + ^{a s - p){ai + p)]' 



(42) 



The respective gain parameters are 



XPL- 



= \fv - \Jv - 1, 



p v = 



1 



l hp 

2^10 feo^p 
1 -Lp~ 
2^10 k w* 



(43) 



To obtain the above decomposition we assumed that: l/w p <C fcoj l/ws,r -C ko, l/w^ r fco, l/kow^, <C 7r 
and Lp/v40 <C u> p . In other words, beam diffraction angles need to be much smaller than the phase matching 
angles, the total walk-off cannot be much larger than the pump waist and the mode arc width has to be 
smaller than n. 

When angular width is relatively small, the modes are just elliptic Hermite-Gaussian functions [25] in spatial 
frequencies. The explicit condition is u> s , r <C kow^, for the signal and w^ r <C kouu"^ for the idler. Then the 
Fourier transform of ( 41 1 yields in the spatial mode functions in the near field 



-ikox 



(44) 
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In particular the i/>oo{x,y) is an elliptical Gaussian beam, traveling along the ideal phase matching cone, 
that is, at the angle arctan(a s ) and passing through the center of the crystal. Its widths are w S)r in the 
radial direction and w v in the angular direction. The mode is optimal in the terms of maximal achievable 
gain, as well as the signal-to-noise ratio. 




FIG. 7: First three modes for the significant walk-off case: tpofi(ks,±, ip s ), ij)™i(ks,±, y>s) and i/>o" 3 (fc s ,x, <p s )- The phase 
e -'k Si „L/2 j g no ^ s j 10wn -phe modes are bent elliptic Gaussian-Hermite beams. When angular width is sufficiently 
small, the mode with the highest amplification -i/>o n o(&s,x, fs) is just an elliptic Gaussian. 



5. NUMERICAL SIMULATIONS 



Above, we have derived the singular value decompositions (|32j) and (37) from the low gain regime scattering 
kernel Si S given in (13). However, along with mathematically well justified approximations we have taken a 



few less rigorous steps, especially (17) and (34). Furthermore, the separation of variables (15) needs to be 

( 33 1 . To confirm the validity of the 



verified, as well as the Taylor series approximation around <p a 
derivation of the characteristic modes, we performed a numerical simulation. 

For each crucial step in the approximation of the scattering kernel Si S , we performed its numerical SVD. 
Technically, we calculated a discretized kernel in a cylindrical coordinate system (18), corrected with the 
proper Jacobian. It was too memory consuming to calculate Si S over the entire rectangular sector of the 
coordinate grid, so we calculated its values only in those regions which have a potential to contribute 
significantly, and we employed sparse arrays. Since Si S ( 13 ) has the pump amplitude ( 16 1 as a factor, 
we took into account only those regions in (fc^Xj &s,x) space for which the values of A p (k s ^± + fcj,j_) are 
significant. In the cylindrical coordinates this happens when 



2.5 



and 



\fs - 'fill < 



2.5 

k w p ' 



(45) 



We chose the constant 2.5 as it covers over 0.999 of the pump intensity. Moreover, we restricted ourselves 
to values close to the ideal phase matching setting fc s ,x, fc-i.x — and fsifi — 0. The extent of the grid 
must be much greater than the zeroth mode, while the mesh must be much smaller than the peak sizes. The 
requirement is especially important for the sine peak, as a too low resolution may spoil its oscillating shape. 
To check if the simulation works properly, we compared results for the same physical data but plotted over 
a twice finer grid or at a twice broader range. 

After verifying the reliability of numerical SVD, we have compared results from different scattering kernel 
Si s approximations for a significant walk-off setting. The results are illustrated in Fig. [8j For this data 
the zeroth singular value £o.o was preserved by the approximation, within a 2% margin of error. However, 
other gain parameters changed visibly with every approximation. The radial shape of the zeroth mode was 
altered only by the sinc-exp approximation (17). Even though the consecutive approximations modified the 
angular width w StV , the modes were still qualitatively Hermite-Gaussians. Ten first characteristic modes 
were verified to be indeed well separable. They could be reproduced by the product of radial and angular 
functions in over 98% of their intensity. 
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FIG. 8: Comparison of i) the gain parameters £o,m and ii-iv) sections though mode functions obtained for four 
different approximations of the scattering kernel Si a : (a) the original with sine (13 1, (b) the kernel with exp (19 1, 
(b') the quadratic form approximation (251, (c) the kernel with the mixed term removed (33 1 and (d) the final 



result (37 1. As (b) and (b') were indistinguishable, they are represented only by (b). Calculations were carried 
out for a /3— BBO crystal with L = 1mm, 9 = 30°, w p = 20/im, X p — 0.4/im, X a = 0.6^im. This corresponds to 
significant walk-off, pL/w p ~ 2.4. The kernels have been calculated on a grid of 61 x 161 points, covering the range 
of [0.6/xm- 1 , 1.2/um" 1 ] x [-§ , §]. 



6. PRACTICAL SIGNAL OPTIMIZATION AND NOISE REDUCTION 



In the previous chapters we worked in the perturbative regime, that is, we assumed that the gain parameters 
are small £„ im < 1. In particular, we found that the gain parameters are proportional to the pump amplitude 
£n,m oc P and that the characteristic modes do not change with the pump amplitude P. However, a useful 
amplifier requires gain parameters much higher than 1. Hence, we need to extrapolate the results. We 
conjuncture that the proportionality relation £ n , m oc P holds also for higher gains, while the mode functions 
remain unchanged. This is strictly true for a 1-D parametric amplifier we considered in Section 2 and was 
numerically verified for a waveguide amplifier pumped with ultrashort pulses |14j . 

In a seeded OPA, two processes occur in parallel: amplification of the seed beam and generation of the 
spontaneous parametric fluorescence. They do not influence one another as long as the saturation effects 
can be neglected. Both of these processes can be described within the framework of the developed model. 
The latter process typically sets the noise level of the amplifier. Since the total amount of fluorescence 
depends only on the pump parameters and not on the seed, it is best to chose the seed beam size which 
provides maximum possible amplification. This is exactly the beam described by the fundamental input 



mode ipofl(ks,-0 gi ven m (41) and ( |44| . In a typical situation, this is an elliptic Gaussian beam of size 
w Str x w v propagating at an angle arctan(a s ) inside a nonlinear crystal. In the interaction with the pump 
light the seed beam is amplified by a factor G = exp(2£ ,o) an d a beam in a mode i>^o(k St x) is produced at 



1G 



the output. 

An elliptical Gaussian beam can be easily converted to a Gaussian beam with the help of astigmatic optics. 
If astigmatic shaping of the input beam is unfeasible, one can use round Gaussian beam of the waist w s — 
y/w s , r Wip- in such case many input modes tp^n 2m(^,i) will be excited, each of them will be amplified by a 
factor exp(2^2n,2m) and a distorted beam will be produced on the output. Typically we may neglect higher 
order modes, since they both have lower gains and are only slightly excited. The most important effect of 
using a round input beam is reduced coupling of the seed light to the fundamental mode. This is described 
by the factor 



r, = --"^ (46) 



which effectively reduces gain to r\G. 



Let us now proceed with the calculation of the intensity of the fluorescence noise emitted per unit time in 
a certain direction. From basic quantum-mechanical considerations we know that the number of photons 
scattered into mode VCm^s.J-) at a frequency of uj s is equal to sinh 2 £ rM „ ~ exp(2£„ jTO ). When we observe the 
output of the amplifier through an aperture, the contributions of distinct modes should be added incoherently. 
Thus, the spectral density of the number of photons is equal to the sum 

(n(w s )) = 2jexp(2f„ )m )T„, m (47) 

n,m 

where T„. m is the transmission of the output mode ^""^(fc^ ,j_) trough the observation aperture. 

To calculate the number of fluorescence photons emitted per unit time d(n) j dt we add contributions from 
all the frequencies 

™ = J du a (n( w> )>, (48) 

where the integral should be performed over the observation bandwidth. With such a formulation, the 
pump intensity and the number of fluorescence photons emitted per unit time d(n)/dt can become slowly 
varying time dependent quantities. However, let us focus on the simplest way of minimizing the total noise 



represented by the integral (48), that is minimizing the noise power at each frequency h(uj s ) independently. 



We consider achieving this goal by spatial filtering of the output of the amplifier with an aperture. 



As one my find in (43) the gain parameters £ n , m form geometric series with respect to both n and m as 
£„ jm cx //V™- For typical data, /i is small and modes with n > can be neglected. On the other hand, 
is often close to 1 and spatial filtering in the angular direction can become necessary. This can be 
accomplished by imaging the crystal onto a suitably oriented slit. One may then achieve high transmission 
for the fundamental mode T^o carrying the useful signal, and significant suppression for the higher order 
mode T m , as plotted in Fig. [9j Results for arbitrary slit width can be checked with a script we provide 

online [55]. The explicit expression for the intensity transmission of tp'^ l u ^ n (k s ,±) mode through the slit of 
width h is given by the integral 

' 2 dyu 2 m (JL) , (49) 



h/2 



where u m (y /w v ) are the angular components of the mode functions (44). For an aperture transmitting 70% 



of the zeroth mode, the signal-to-noise ratio may be elevated up to 3.3 times. 

Another approach to reducing fluorescence contained in high order modes is reducing angular gain decrement 
Our simplified model predicts that it would approach zero when the following relation holds between 
the crystal length and the pump beam waist 

(50) 
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FIG. 9: Plot of intensity transmittance To, m of the first 15 Hermite-Gaussian modes trough a slit set to transmit 
70% (squares) or 90% (dots) of the zeroth mode. 



Naturally the above relation should be considered approximate, because it relies on the numerical factor 



used in sine approximation (17). Fortunately, even a rough fit to the (50 1 should suffice 



CONCLUSION 



We have developed a simple analytical model of a noncollinear parametric amplifier pumped with a focused 
monochromatic beam and utilizing type I phase matching. We found an approximate Bloch-Messiah re- 
duction for a low-gain parametric amplifier for an arbitrary pair of signal and idler frequencies. The final 
result was expressed in a closed form in two special cases: zero walk-off of the pump beam or total walk-off 
of the order of the beam size. Characteristic modes of the first setting are optical vortices. The latter 
case corresponds to the typical experimental situation and we find that the characteristic modes are elliptic 
Hermite-Gaussian beams. We checked the validity of the approximations assumed during the derivation by 
comparing the analytical result to the numerical calculations for a typical range of parameters. Finally, we 
scaled the results of the reduction to a high gain regime using the results of our previous work |14| and we 
discussed how to calculate the gain and fluorescence intensity. 

In particular, we have calculated the fundamental mode of the amplifier, which has turned out to be elliptic 
Gaussian. Seeding the amplifier with a beam matching this mode yields the highest possible amplification. 
The output of the parametric amplifier contains both the amplified seed beam and optical noise due to para- 
metric fluorescence. The Bloch-Messiah reduction allows us to directly calculate the amount of parametric 
fluorescence emitted by the amplifier. Further analysis shows that it can be reduced by either suitable spatial 
filtering of high-order modes or adjusting the pump beam diameter so that that those modes are suppressed. 

The results of our approximated model will never be as accurate as numerical simulations based on tracing 
the evolution of quasi-probability distributions |18) . but they can be instantly calculated for various con- 
figurations of parametric amplifiers using the Mathematica 6 script we provide online [26] . and thus may 
provide valuable insight in the applications. 

Our results may be also used as a starting point for numerical code modeling OPCPA with pump depletion. 
This would be accomplished by dividing the amplifier longitudinally into two parts: the front one, in which 
a gain of 10 3 -10 6 is reached without depleting the pump, and the rear one, where saturation occurs. In the 
first part, quantum phenomena play a role, while in the second part spontaneous fluorescence can be treated 
as classical noise. The front part can be treated with our model. This way the optimal input seed beam 
shape and spatially resolved intensity of the parametric fluorescence can be obtained. The latter can serve 
as an initial condition for the numerical code solving the classical evolution of the fields in the rear part of 
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the amplifier. 

The model developed in this paper may be also helpful for optimizing photon pair sources, since it represents 
an alternative approach to the problem of finding spatial modes of the parametric fluorescence and optimal 
fiber coupling |16) . 
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